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Abstract: In backgrounds with compact dimensions there may exist several phases of 
black objects including the black-hole and the black-string. The phase transition between 
them raises puzzles and touches fundamental issues such as topology change, uniqueness 
and Cosmic Censorship. No analytic solution is known for the black hole, and moreover, 
one can expect approximate solutions only for very small black holes, while the phase 
transition physics happens when the black hole is large. Hence we turn to numerical 
solutions. Here some theoretical background to the numerical analysis is given, while the 
results will appear in a forthcoming paper. Goals for a numerical analysis are set. The 
scalar charge and tension along the compact dimension are defined and used as improved 
order parameters which put both the black hole and the black string at finite values on 
the phase diagram. Predictions for small black holes are presented. The differential and 
the integrated forms of the first law are derived, and the latter (Smarr's formula) can 
be used to estimate the "overall numerical error". Field asymptotics and expressions for 
physical quantities in terms of the numerical ones are supplied. Techniques include "method 
of equivalent charges", free energy, dimensional reduction, and analytic perturbation for 
small black holes. 
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1. Introduction 

In the presence of extra compact dimensions, there may exist several phases of massive 
solutions of General Relativity, depending on the relative size of the object and the relevant 
length scales in the compact dimensions. For concreteness, we consider a background 
with a single compact dimension - M'^"^'^ x S^. We denote the coordinate along the 
compact dimension by z and the period by L. The problem is characterized by a single 
dimensionless parameter, for instance^ the dimensionless mass ix = GnM/L'^-^ where Gn 
is the d dimensional Newton constant and M is the mass (measured at infinity). 
The relevant phases are 
^ Later we will define an alternative parameter x. 
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• The uniform black string. 

A string solution is one which has an S'^~^ x horizon topology. The uniform string 
is ^-independent and it is described by the d — 1 Schwarzschild - Tangherlini metric 
[Q] with the addition of +dz'^ - a spectating z coordinate. This solution was shown to 
be unstable to gravitational perturbation below a certain critical by Gregory and 
Laflamme (GL) [|, |. 

• The non-uniform string. 

These are solutions with the string horizon topology, which are no longer z-independent 
The marginally unstable GL - mode implies that such a branch of solutions is ema- 
nating from the GL critical point (for instance due to Morse theory [Q]), and a priori 
there could be other branches as well. Gubser determined that the phase transi- 
tion at the GL point is first order by an analytic perturbation analysis |5|. Morse 
theory implies that the emanating branch would be unstable, at least at the be- 
ginning Q. Finally these solutions were found numerically by Wiseman in 6d Q, 
who managed to formulate axially-symmetric gravitostatics (namely, essentially 2d) 
in a "relaxation" form (a procedure familiar from electrostatics) while presenting the 
constraints through "Cauchy-Riemann - like" relations. 

• The black hole. 

These are solutions with S*^"^ horizon topology. For small (dimensionless) mass we 
expect these solutions to resemble a d dimensional black hole near the horizon. No 
analytic solution is known, and we consider our numerical solutions in 5d as strong 
evidence for their existence (indications appeared already in Q). 

• The stable non-uniform string. 

Horowitz and Maeda |^ postulated the existence of an additional stable phase, to 
serve as the end-point for the GL decay. This phase should be distinguished from the 
unstable non-uniform string which is presumably a consequence of Gubser's analysis 
(or from strings which are too massive to serve as the end-point for decay - clearly 
the stability of the solution is critical for the physical interpretation). Horowitz 
and Maeda argued, based on theorems for the increasing area of the event horizon, 
that horizon pinching is impossible. By now there is mounting evidence against the 
existence of this phase: in Q it was shown that a continuous topology change (as a 
function of a parameter rather than time) is actually possible, moreover Morse theory 
implies an obstacle, namely that the addition of this phase must be accompanied by 
some other one, in ^ it was shown numerically that the whole branch of static non- 
uniform strings originating from the GL point is too massive to serve as an end-point, 
finally a full time-dependent numerical simulation had to stop at -Rj^jj^Z-Rmax ~ 1/13 
(due to grid stretching) without finding clear evidence for stabilization towards this 
phase |]lO| ^ (i?max, -Rmin refer to the z-dependent radius of the horizon) . Additional 

^Note however that the authors of this paper did not interpret their results either as supporting or as 
countering the conjecture. 
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circumstantial evidence is supplied by Gubser's frustrated hope to find a second order 
transition H] and by the failed attempt to find the non-uniform string analytically 
by looking for 5d algebraically - special solutions [p!l| ]. 

Our interpretation is that the argument of |p fails at a deep and interesting level - 
the horizon theorems being used there require the singularity not to leave the horizon, 
while it seems to us plausible that this system constitutes an example for violation 
of Cosmic Censorship. 



In addition to the research described above we would like to mention: an ansatz that 
reduces the number of unknown functions due to Harmark and Obers speculations for 
a generalization of uniqueness to d > 4 [13|, some phenomenological implications for huge 
explosions accompanying the transition [0] , numerical evidence for the ideas of Q includ- 
ing the approach to a cone |15, 16], a time-symmetric initial data study |^], a numerical 
study of black holes on a brane which have a similar geometry |17] a recently discovered, 
and a related instability in highly rotating black holes [18|. Other related research included 
a relation between the GL transition and gauge theories through Matrix M-Theory [19|, a 
prediction for a non-uniform near-extremal black string [^], a useful review of numerical 
relativity [^], a relation between the GL instability and the negative heat capacity of black 
holes 1 22 1, a discussion of implication of the GL instability to more general spacetimes with 
horizons [23|, a discussion of the explosion expected during these phase transition |24], a 
recent analysis of 4d metrics with a compact dimension but different asymptotic conditions 

I |29|1 for 



and a study of non- uniform strings in the two-brane system [26|. See also |2^, |2^, 
a short and non-representative list of papers on black holes and large extra dimensions/ 
brane worlds/ accelerator prospects. 

Here, motivated by Q we take the route of approaching the phase transition by in- 
creasing the mass of the black hole, rather than through the GL instability. Since no 
analytic solution is known for the black hole we turn to numerical solutions. Moreover, 
even though one can expect approximate analytic solutions to exist for very small black 
holes, the phase transition physics happens when the black hole is of the same size as the 
compact dimension, and so a numerical simulation is essential. It is important to list the 
goals of such a numerical study, even if currently we cannot reach all of them 



• Establish the existence of the black hole solution. 



• Determine the maximal (dimensionless) mass for a black hole (it is expected not to 
fit in the compact circle beyond a certain size Q). Identify the direct reason for that. 

• Establish another failure of generalized uniqueness, by demonstrating a black hole 
solution with the same mass as some (stable) black string. 

• Test the description of approach to topology change proposed in Q , namely that the 
geometry in the region between the two poles approaches a cone and that at least for 
5 < d < 10 the geometry destabilizes (perturbatively) before reaching transition. 
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• To determine whether the black string decay violates Cosmic Censorship is one of 
the most interesting issues pertaining to the system. However, it is not clear if the 
study of black holes would contribute directly. 

This paper contains theoretical background and results related to the black hole so- 
lutions. After describing our ansatz for the metric in section |2| we turn to an analysis of 
the asymptotics and charges at infinity. We start in section ^ (inspired by ||3^, 31]) with 



what we term "a method of equivalent sources" where one pretends that gravity is weak 
throughout and the source is large and smeared. We identify the numerically measured 
asymptotics both in harmonic and in conformal gauge and find their relation with the 
physical charges, the mass and the tension. We promote a physical picture of the tension 
being a tension of an imaginary string, which acts to decrease the size of the extra dimen- 
sion counteracting the effect of the mass which "wants to make space" for itself. We also 
stress the importance of the tension as a natural coordinate for the system's phase diagram 
instead of the previously used A = ^[(i?max/-Rniin) ~ ^]) "-"^^ ^^^^ advantage of 

putting both the black hole and the black string on the same phase diagram, being always 
finite. This section ends by making some predictions for small black holes. 

In section ^ we use a rigorous technique based on the (Gibbons-Hawking) free energy 
to revisit some of the issues of the previous section and to derive the first law. Together 
one arrives at a form of the first law, dm = T dS + t dL which is completely analogous to 
the familiar dE = T dS — P dV (for gas) (see also ) . This (differential) form of the first 
law is then readily combined with scale invariance to yield the integrated form of the first 
law (also known as "Smarr's formula"). The latter is an important test for numerics since 
it ties quantities at the horizon with those at infinity and relies on the satisfaction of the 
equations of motion everywhere (since the derivation uses integration by parts) and so it 
gives a measure for the "overall numerical error" . 

Next, in section ^ we use the fact that in the asymptotic region the z dependence is 
lost to perform a dimensional reduction. The mass defined in section ^ is seen to agree 
with the standard definition of mass in the lower dimension. The asymptotic constant for 
the decay of the size of the extra dimension is interpreted in the lower dimension as a scalar 
charge. The scalar charge, like the tension is a natural parameter for the system's phase 
diagram and it is especially useful for the Gregory-Laflamme transition since it vanishes 
for the string and is non-zero for all other known solutions. It shares another interesting 



property with the tension: the tension, much like the mass, is always positive |33, ^ and 
it seems that so is the scalar charge ~ it would be interesting to prove/ refute this. 

In section ^ we turn to the horizon, identify the measurable quantities, and give pre- 
dictions for small black holes. We conclude by summarizing our main results in section |^. 
The appendices include details of a coordinate transformation, yet another confirmation of 
the expression for the mass (this time using the Hawking-Horowitz expression), the details 
of action evaluation in the harmonic gauge and a different proof of Smarr's formula. 

Very recently the paper appeared which overlaps with some of these theoretical 
considerations. Our 5d numerical results will follow in a sequel paper Since the 
appearance of the first version of this paper, another interesting paper appeared |36|. 
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2. Ansatz for metric 



The system under study, a static (no angular momentum) black object in M''"^'^ x S^, is 
characterized by three dimensionful constants: L the size of the extra dimension such that 
z and z + L are identified, M the mass of the system measured at infinity of M*^"^, and Gn 
the Newton constant in d dimensions. From these a single dimensionless parameter can be 
constructed 

^i = {GnM)/l^-'\ (2.1) 

while the {d — 1) dimensional effective Newton constant is given by Gd-i = Gn/L. 

The isometrics of these solutions are 0{d — 2) x U{1) x (^2)^) where the 0{d — 2) 
comes from the spherical symmetry in M'^^^ which we assume and the ?7(1) comes from 
time independence.^ The discrete factors represents reflections with respect to time and 
the z coordinate t —t, z — > —z. The most general metric with these isometrics is 

ds^ = _e2^dt2 ^ds^^^^^ + e2^d02, (2.2) 

which is a general metric in the (r, z) plane together with two functions on the plane 
A = A{r, z), C = C{r, z). The horizon is a line determined by e^^ = and dO^ is the line 
element for S'^"'^. 

There are several alternatives to fix the gauge and choose coordinates in the (r, z) 
plane, though for most of the current theoretical analysis this does not matter. For the 
numerical convenience we choose the ansatz 

ds'^ = -A^dt"^ + e'^^{dr^ + dz'^)+r^e'^^ dn'^ (2.3) 

where the conformal gauge is chosen for the 2d metric and the new variables A, G are 
defined in terms of the variables in eq. ( |2.2[) 

A:= 

G :=&- log(r) (2.4) 

so that A is regular at the horizon (^|}^or ~ ^) ^"^^ such that C ^ at infinity. 

This ansatz still contains a gauge freedom for conformal transformations which we 
finally fix by setting the domain of definition for the three functions to be {(r, z) : \z\ < 
L, + > where we define L := L/2 for numerical convenience. It was shown in 



1 17] that transforming to this domain is always possible by writing down elliptic equations 
for the coordinate transformation. Thus each domain is characterized by a dimensionless 
constant (the conformal invariant) 

X := 2ph/L. (2.5) 

The normalization was chosen so that 2; = 1 is the maximum possible value. Since the 
problem has only one dimensionless parameter all quantities will be a function of x (and 
in this sense it replaces the dimensionless mass ^). 

We still need to state the boundary conditions. There are four boundaries 



^In the Lorentzian solutions it is the non-compact version of (7(1). 
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• r — oo infinity. 

The metric is asymptotically flat 

Qt^v ~ r]^^ := diag(-l, 1, . . . , 1) <^ 

= 1, S = C = (2.6) 

• Horizon - a regular horizon. 

• r = axis - regularity on the axis. 

• \z\ = L/2 ~ periodicity (and reflection symmetry at z = 0). 

By comparing with the Schwarzschild metric one can flnd the relation p = p{ps) where 
ps is the Schwarzschild coordinate (see appendix ^ 



P_ 

Ph 



Pa \ Po 



2 

d-3 



(2.7) 



where ps = po is the location of the horizon in Schwarzschild coordinates. Asymptotic 
flatness chooses a natural radial coordinate at infinity, hence we require that asymptotically 
p ^ Ps, and we find the relation between ph and po 



2 



Ph = po/2—K (2.8) 
3. Charges and the method of equivalent sources 

We start by analyzing the system at r — >■ oo. In this limit all the metric functions become 
independent of z (since the mass of the ^;-dependent modes is proportional to 1/L all z 
dependence drops like exp(— r/L) as one gets away from the black hole). Moreover, at 
infinity we are in the weak field limit, namely the Newtonian limit. 

It is possible to derive the properties of the asymptotics by pretending that the black 
hole source is smeared over a much larger area such that the gravitational field is everywhere 
weak, while the energy-momentum tensor is non-vanishing (recall that there are no sources 
in the sought-for black hole solutions). This method is inspired by the analysis of pO, 31 1, 



and should be valid since far away from the source it should not matter whether the source 
is a black hole or a low density star. We term this approach "a method of equivalent 
sources". After we derive the expression for the charges in terms of the constants of 
asymptotics within this method, we return in the next section and prove it rigorously 
using the free energy and the first law. 
One defines 

hf_ii/ '■= h^i, '^h^^T^fii/ (3-1) 
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where r/^,y is the flat space metric, which is used for raising and lowering indices, and h <^ 1 
is the perturbation. Choosing the harmonic (or "Lorentz" ) gauge 

9^ V = 0' (3-2) 

the linearized field equations become 

Ah^u = -IGnGNT^iuix) (3.3) 

where T^y[x) is the stress energy tensor. The solution may be written in a multipole expan- 
sion.^ The leading term at infinity comes from the monopoles (dipoles, including angular 
momentum effects vanish here due to the reflection symmetries, so the next correction will 
come from quadruples), and in our geometry^ it is 

Ku = 2 kd r£ (3.4) 

where 

7^?. := / d^-^xT^, (3.5) 
are "total mass-energy charges" at t = const, and 

■■= p- (3.6) 

■= V/2 + i)=^(d7^- ^'-^^ 

firf-i is the area of the unit S'^~^ sphere, and for odd d we define {d/2y. := . . . |. 

From now till the end of this section we will use "asymptotics adapted units" such that 
Gd-i = Gn/L = 1. 

Finally one inverts the transformation to retrieve h 

h^u = — ^ _ 2 ^oH^iv — 

T^u-:r^V,u], (3.8) 



{d - 4) \ f"' d-2 



where := T^^??"^. 



3.1 Mass and tension 



The symmetries of the problem restrict the form of the energy- momentum charges T^: 
time reflection implies = = 0. Similarly, z axis reflection implies = 0. Spherical 
symmetry implies that = 6ij. Here i,j are any index other than = t,z. 



^where the constant of proportionality is read as usual from A(l/r'*^*) = {d — 4) S{x). 

^In order to get the solution in uncompactified spacetime replace (rf — 4) Lr'^^'^ (d — S) fld~2 ' 
in the following formulae. 
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It is important that not all three energy- momentum charges T^, T^z are indepen- 
dent, as a result of GR's constraint equations. The relation we are about to derive is a 
consequence of the constraint Grr = (where G is the Einstein tensor), but here we show 
how to get the same result from conservation of energy-momentum 

S'^T^. = (3.9) 

where we used = 9^ + 0{h). In particular 

5'^ Trr = Q^ Trr{r) = const (3.10) 



but 



and hence 



We denote 



Trr{r ^ L) = ^ Trr{r) = (3.11) 
T^^ = 0. (3.12) 
m := , (3.13) 



the mass of the black hole. Tzz is pressure and since the conservation law tells us it is 
constant along z: dz Tzz = it makes sense to define 



r := 



-TQ/L (3.14) 



to be the tension along the periodic direction. This is the tension measured by an asymp- 
totic observer, which interprets it as the tension of an imaginary string stretched along 
the compact circle. The choice of sign reflects the fact that the tension is always positive 
1 33 , 34 1 . Altogether we write the energy- momentum charges of the black hole as 

r« =diag(m,-TL,0,...,0) (3.15) 

3.2 From harmonic to conformal coordinates 

The most general ansatz for the asymptotic region in harmonic coordinates ( |3.2| ) is 

ds^ = -e-2^« df + e"^" dz^ + e2^« {dr^ + dnj_,) (3.16) 

where Ah, Bh, Gh depend only on the harmonic radial coordinate rn, and the sign for Ah 
was defined for later convenience. 

Define the constants an, bn, ch (with dimension length'^"^) from the asymptotic form 
of the functions 



Bh = -£j + 0{ 



Ch = ^,+0{^,). (3.17) 
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The harmonic gauge condition (|3.2| ) sets some relations between these constants. Ac- 
tually, only d^hrr is not identically zero and gives 



an + bn + (d - 4.) ch = 



(3.18) 



Eqs. (U, |rT|) tell us that 



an 
bH 

CH 



kd 



{d - 4) (d - 2) 



d-3 -1 
1 -(d-3) 
1 1 



m 
tL 



(3.19) 



Note that owing to this relation the constraint ( 3.18 ) is satisfied automatically, which can 
be traced to the conservation of energy- momentum ( 3.1C| ). 

Comparing with the conformal ansatz (2.3), where this time the functions vl, B, C 
depend only on r, one gets 



A = e-^- 
B = Bh 



dr 



drn 



C = Ch- \og{r/rH) 



(3.20) 
(3.21) 
(3.22) 
(3.23) 



The equation for r from (|3.22|, 3.21|) 



dr 



drn 



I drn 



(3.24) 



can be integrated to give 



for d 



(3.25) 



where ro, ruo are constants of integration which are needed in the numerics in order to 
keep the location of the horizon fixed. Note that the case d = 5 is somewhat special and 
needs to be treated separately. 

We are now ready to find the constants in the conformal gauge asymptotics. For ^4, B 
the work is done already due to ( p. 19 - 3. 21 ). We define a, b by 



I- A 
B 



a ( ^ \ 
1 



^d-A 



and we get 



a 






b 




bH 



kd 



{d - 4) (d - 2) 



d-3 -1 
1 -{d-3) 





m 




tL 



(3.26) 
(3.27) 
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which can be inverted to give 





m 


_ j_ 


d-3 -1 




a 






tL 


kd 


1 -(d-3) 




b 





(3.28) 



(reminder: is defined in (p.6|)). 

The last equation is one of our main results - it tells us the physical charges in terms 
of the numerical asymptotics. 

We now turn to the function C(r) which we treat separately for d = 5 and for d > 5. 
For d = 5 an expansion of ( p. 23] ) using the change of variables ( p.25| ) gives 

CH + hrL log(rj^o) ^tL 



C 



log(r) + 0(i^; 



/yt /yt 

The leading term at infinity is log(r)/r and so we may define C5 by 

C = C5 h... , 



and using (|3.2g|, |3.28D we get 



C5 = 6h5 - = 2 6 - a. 



(3.29) 

(3.30) 
(3.31) 



This constraint is the 5d manifestation of ( 3.1^ ). Note that the next term 0(l/r) is hard 
to distinguish from the leading term and it contains an arbitrary constant rno- 
For d > 5 (|]2|) yields 



kdTL 



, fj ru {d-4:){d-5)r 

which together with the inversion of ( |3.25| ) 

kdTL 



d-4 
H 



+ 0{- 



d-5 



)) 



H 



ro + 



+ 0{- 



.d-i' 



(3.32) 



(3.33) 



{d - 4) (d - 5) r'^-s 

tells us that the leading term in C(r) is — while the charge r enters only at order 0(^3^4 
and is mixed with tq. 

3.3 Small black holes 



Myers |31| argues that for small black holes (x <C 1) one should take zero tension, r = 0, 



namely the equivalent source is dust. In this case we use (3.27) to compute the asymptotics 

(3.34) 





a 


r^j /m 


d-3 






b 




1 





where the ~ sign denotes that the ratio of the two sides approaches 1 as x — > 0. 
Moreover, in this limit the mass is determined by the horizon size 



-1^^+^ m ~ if)'^-^ = 4 (f ) 



d-3 



{d-2) Ld-3 



(3.35) 



This can be seen by using the method of equivalent sources to compute the mass of a 



Schwarzschild black hole |3C] and then noting that after compactification the mass does 
not change in the weak gravity region since the Newtonian potential satisfies Gauss' law. 
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4. Charges, free energy and the first law 



4.1 Charges from free energy 

Let us derive the form of the first law (conservation of energy) using rigorous free energy 
techniques. We will see that the results agree with those of the previous section. 

The (Gibbons-Hawking [^]) free energy action integral is defined for the Euclidean 
continuation of static metrics and is given by 

^ = = TTT^ / '^^d.^d + / dVd^i [K - K^] (4.1) 

IQ-kGn J 8itGn Ja 

where /? = 1/T is the inverse temperature (the period of the Euclidean time), F is the free 
energy, Ra is the d dimensional Ricci scalar, dVd, dVd-i are the d, {d — 1) invariant volume 
elements, the second term is integrated over the boundary of the manifold d, K is the trace 
of the second fundamental form, while is the same quantity for a reference geometry 
with the same boundary. 

Consider the action to be a function of the boundary conditions 



I = I{P,L) 



(4.2) 



evaluated after extremization over metrics. In order to get the first law we need an explicit 
expression for the action, find the total differential dl = {dl /d(3) d[3 + {dl /dL) dL and then 
use thermodynamic identities to transform it into the standard form. 

The action in terms of the harmonic fields ( |3.16| ) reads (see appendix ^ for the deriva- 
tion) 



2k^ Gn 



where we use a short-hand notation Xi 





K 



[Ah-, Bh, Ch] and the matrix Kij is 
-1 



-(d-3) 
-1 {d-3) 

(d-3) (d-3) {d-3){d-4) 



(4.3) 



(4.4) 



Moreover, there is no boundary term. 

Having the expression for the action we wish to compute its differential. By a general 
theorem in mechanics ^ it is given by the conjugate momentum at the boundary 

dC , 



dxil 



drXi 



\r„=R 



(4.5) 



where C is the Lagrangian density. 

Using the constants of asymptotics ( p. 17 ) we may compute 













I 




dcH _ 





JoL_ 
kdGN 



id-4:)K 



an 
bn 

CH 



(4.6) 



®See for example |^^: 5/ = Jgj + J^^i^ — ) '^9'^^ since the second term vanishes as we 
are evaluating on solutions of the equations of motion, only the first term contributes. 



- 11 - 



where ( [4 .51) was evaluated at the i? — > oo hmit. 

is precisely proportional to the constraint (|3.18| ) (derived from Grr = 0) and thus 
vanishes identically. This is to be expected since a shift in C is equivalent to a shift of the 
boundary in ru without changing the asymptotic sizes /?, L, and since we are evaluating 
/ in the R ^ oo limit, the action should not change. Now we may use the constraint 
(3.18) to eliminate ch from the expressions ( [4.6D and we may drop the H subscript from 
a, b according to ( |3.27| ) and get 



(4.7) 







d-3 -1 




a 


dB„ _ 


kd Gn 


-1 d-3 




b 



Finally the relations dp = —PodAn, dL = L^dBn valid asymptotically where Ah 
Br = yield the differential of the action 



dl = -{{d-3)a-b)dl3-i{a-{d- 3) b) dL 
L L 



(4.8) 



Now let us compare with standard thermodynamics. In analogy with gas thermody- 
namics where one has 

dE = TdS - PdV (4.9) 



the first law must take the following form 



dm = T dS + T dL 



(4.10) 



This can be considered the proper definition of the mass m and the tension r. We will 
now see that this definition coincides with the definition in terms of asymptotics ( ^.28 ) and 
hence all is consistent. The free energy is related to m by -F = m — TS* and the action is 
/ = —(3 F. Putting the last two relations together one finds 



dl 



-m dp — (3 T dL 



(4.11) 



Comparing with (4^) we re-obtain the expressions for the mass and tension ( 3.2^ ) as 
anticipated above. Alternatively, if one takes (|3.28| ) to be the definition of mass and 
tension, then (^^) proves the form of the first law ( [4.10| ) and suggest the analogy with gas 
thermodynamics (4.9). 



4.2 Smarr's formula 

In this section we describe Smarr's formula, also known as the integrated form of the first 
law, for the geometry under study. It is a relation between thermodynamic quantities at 
the horizon with those at infinity, relying on the generalized Stokes formula and the validity 
of the equations of motion in the interior, and as such it estimates the "overall numerical 
error" in our numerical implementation |^. 

In general Smarr's law can be gotten from the first law ( [4.10D after taking account of 
the scaling dimensions. Performing a scaling transformation L — > (1 + e)L, recalling the 
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relations T = h/{2tt), S = A/{AGm), and that the area has dimension d — 2 while the 
mass has d — 3, we get by we expanding ( 4.10I ) to first order in e 



{d-3)m = id-2)^ + TL 



(4.12) 



In order to put ( [4.12 ) into a "numerically adapted form", that is one that uses only 
the numerical quantities we use ( |3.28D to get 



1 



8 TT GjV 



-Ak 



{d - 4:)L 
kd Gn 



a, 



(4.13) 



where on LHS we have horizon quantities and on the RHS asymptotic ones. Finally the 
LHS can be written also as T S. 

5. Dimensional reduction and the scalar charge 

At infinity, since all z dependence is lost, it is natural to look at the system from a (d — 1) 
point of view. We will see that this point of view offers a natural interpretation for the 
mass formula (|3|) - the mass can be read from the metric in the usual way. In addition, 
this perspective suggests the importance of the scalar charge to be discussed next. 
The d dimensional metric can be written as 



ds^ 



-g^^-^Ux'^dx' + e^Uz^ 



(5.1) 



where a, h run over the d—\ dimensions. Thus after reduction we get a metric and a scalar 
field (f) (for a general dimensional reduction one gets a vector as well, but it vanishes in our 
static ansatz). 

Computing the action after reduction we have (|C.1|) 



lid) 



+ ... 



(5.2) 



where the R is the Ricci scalar and . . . could include other terms such as (dcp)'^, A (</>). One 
notices that the gravitational action behaves as e'^ \/—g^'^~^^ R^'^~^^ and is not canonically 
normalized. In order to mend that a Weyl rescaling is required 



The changes are ( C.3 ) 



-Jd-i) 



(5.3) 



-2w 



+ ... 



-9 



(d-l) 



-9 



(d-l) 



(5.4) 



where the . . . denote this time terms with {dw)"^ and A(w). We see that in order to get 
the Einstein action we should choose 



w 



d-3' 



(5.5) 
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Let us look at /iqq 



^ d-3 
But this should equal {po/r)'^~'^ hence 



^ («-T^^) (5-6) 



2 (a - 6) = (po)"^' = m (5.7) 



where the second equality comes from ( p. 35 ) and altogether we reproduced ( p. 281 ) 



Let us reinterpret the asymptotic constant b in light of the dimensional reduction. 
Since by definition Bh = (f), b describes now the fall-off of the scalar field 4> 

(5.8) 



Such a quantity is usually called the scalar charge (of (p). In 1 40 1 it was shown that the 
scalar charge appears in the first law as {dm/d(j)^^) = Gij S-', where T,^ are scalar charges 
defined in analogy with b in ( ^.8D , the indices i,j allow for several scalar fields with a non- 
trivial metric read from the kinetic term Gij{(p) 50* d(j)^ . Indeed after the Weyl rescaling b 
replaces r in the expression for dl (as in ( |4.11| ), while at the same time we must use the 
transformed f3. 



We can now interpret the equation for b ( p.27 ) as telling us that the mass increases b 



and tries to open up the extra dimension, while the tension counteracts. For the uniform 
string these two effects exactly cancel each other. Hence 6 is a useful order parameter to 
describe the departure from the uniform string. Moreover, it has the advantage of being 
finite also for the black hole and hence one can naturally put all the phases on the same 
phase diagram as advocated in Q. We notice that b is always positive as far as we know. 
It would be interesting to understand/ refute this. 



6. At the horizon 



It turns out that all the measurable quantities which we define reside either at infinity or 
at the horizon. After discussing the former in the previous sections, we now turn to the 
latter. In the vicinity of the horizon it is convenient to transform the (r, z) plane into polar 
(p, x) coordinates through z = p cos(x), r = p sin(x). 



6.1 Definition of measurables 

There are two thermodynamic quantities defined at the horizon 
• The area, A. 
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In the conformal gauge ([2. 3D it is given by 

/"TT 

-4 = / ph dx ■ ^d~3 (Ph e'^ sm{x)) 
Jo 



^^-3 r e""^^'-'^ ^ sin(x) dx (6.1) 

Jo 



Surface gravity, k. 



Since the geometry is static we may use the definition of kappa by taking the Eu- 
cHdean continuation of the metric and finding the period /3 of Euchdean time such 
that there is no conical singularity at the horizon. Then k = 2 71/(3. In conformal 
coordinates it is given by 

K = e-'^dpA (6.2) 

The same expression is gotten from the equivalent definition = — {l/2){d^£/^){dy£^^) 
where is the Killing field and the expression is evaluated on the horizon (see for 
example ||4l[|). 



In addition we monitor two quantities that measure the geometry of the horizon 

• The eccentricity, e. 

The spherical symmetry of a small black hole, SO{d — 1), is broken to SO{d — 2), as 
it grows. In the {r, z} plane we define the eccentricity 

1 

A|| := nd-^e^^-""^^ p^-\ at z = 0, 

Ai_ := / Ph dx ■ (e^ Ph sin(x)) '^"^ = 

^0 

= n,_4 pf r e^'-^'-'^'^ sin(x)'^-^ dx (6.3) 
Jo 

A^^,A± are sections of the horizon: is the area of the equatorial sphere at z = 0, 
while A± is the area of the axial (or "polar") sphere at x = 0, vr. 

• The polar distance, ^poles' 

This is the proper distance between the "north" and "south" poles along the r = 
axis and it is given by 

»L/2 

dz 



■^poles •= ^ / dze'^, at r = 0. (6.4) 



^This definition diflters from the standard definition of an eUipse's eccentricity. It is analogous to a/b—1 = 
(1 — e^)~^^'^ ~ 1 + 0.5 e^, where e is the conventionally defined eccentricity. 
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6.2 Predictions for small black holes 

Out of the two thermodynamic quantities A, k one can form a dimensionless quantity 

(6.5) 

For small black holes we may take the Schwarzschild metric as an approximation, 
where A = Po~'^^ = {d — 3)/{2 po), and hence 



d-3 



(6.6) 



More generally we may wish to expand A^'^^ = A^'^\x) in a Taylor series 



AW(x) = ^A 

j=0 



(6.7) 



Note that A^-^"^ depends both on j and on the dimension d which is suppressed in this 



notation. 



It can be shown |42] that the next non- vanishing term is A^J^}^ and it is given by 



(d-2)C(d-3) 



L 



d~3 



(6.8) 



where C is Riemann's zeta function, and we used eq. (|2.8| ). 
For the eccentricity the prediction is l43] 



€ = Cf X 



d-l 



while the constant of proportionality is currently computed only in 5d 



c,((i = 5) = fC(4) = i%=. 2.89 



(6.9) 



(6.10) 



All 5d constants stated in this subsection are confirmed numerically 



7. Summary 



We summarize our main results: 



Asymptotics. 

There are two constants of asymptotics. In the harmonic ansatz ( p.l6| ) there are 3 of 
them defined by ( |3.17| ) and they satisfy the constraint ( ^.181) . In the conformal ansatz 
( |2.3D the constants of asymptotics, a, b, are defined in (|3.26| ), related to the harmonic 
constants in ( 3.27| ), and expressed in terms of the system's mass and tension (around 
the z coordinate) in ( |3.28 ). We stress the physical interpretation of the tension. 
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5d is somewhat special and we can furthermore define the asymptotic quantity C5 in 
( |3.3C1| ), and express it in terms of a, b in ( |3.31 ). 



For small black holes we have a prediction for all quantities in terms of x in (3.34 



Thermodynamics. 



We derive the (differential form of the) first law in ( 4.10 ) and the integrated form in 
&■ 
• Horizon region. 

We make predictions for small black holes for two dimensionless quantities: Ak'^~'^ 
( |6.8D and the horizon eccentricity ( [6.101 ). 
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A. Schwarzschild in conformal coordinates 

Let us transform the Schwarzschild metric into conformal coordinates and find the relation 
between the radii po and ph- The Schwarzschild solution is 

ds^ = -f df + dpi + pI dnl^ 

/ \ d—3 

f = l-0^) (A.1) 



.PS, 

One considers a new radial variable p = p{ps) such that 

rUpi = e^''dp' 

Pi = e'''p' (A.2) 

in terms of which the metric becomes 

ds^ = -fde + e^^{dp^ + p2 dnd-2) = 

= -fdt^ + e^^idp^ + p^ {dx^ + sm\x) (A.3) 
^From ( [A.2| ) one extracts the equation for p 

dp ^ ^-y,d_ps_ 
P PS 

The integral on the right can be solved by a change of variables cosh(/3) = {ps/ Po)^^^'^^^'^ 
and we get 

2 

p/ph = exp[-^— ^arccosh((p5//0o)^'^"^^/^)] = 



Po \ Po 



2 

d-3 



(A.5) 
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The metric in the new conformal (or isotropic) coordinates takes the following simple 
form 

ds' = - (i^)' dt' + (1 + [dp' + p' dnU 

B. Hawking-Horowitz mass 

Let us compute the mass directly from the Hawking-Horowitz formula as an additional "in- 



dependent" derivation. The Hawking-Horowitz formula |43| for a static vacuum solution 



1 

ruHH = -0 (B.l) 

O TT Lr]\f 

where S^'^-^) ^j^g _ 2) 

area of the boundary of a t = const slice, "S^"' is the same 
quantity for a reference geometry, and dn is a derivative with respect to orthogonal motion 
of this boundary. 

We evaluate ( |B.l| ) in the harmonic ansatz ( |3.16| ) for a boundary at rn = R and take 
the R ^ 00 limit. The first term is 

= e-^" drAiLe''")nd-3irHe^")^-^] \r^=R= 

^ L n^.s [i^-^)id- 3) R^-' + {bu + (d - 3) ch)] = 

= Lnd-^[{d-2,)R''-^ + bH]. (B.2) 

The reference t = const geometry is flat: ds^ = dz^ + dr"^ + (irirf_3 with period 
= L e^" and r = R e^" to match with the boundary. The second term is 

~ L Orf_3 (d - 3) + (6h + {d-4) ch)] = 

= LQrf_3(d-3)[i?'^-4 + aH]. (B.3) 
where in the last equality we used the relation ( ^.18 ). 



^This mass definition coincides with the ADM mass when both are applicable. 
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Combining ( B.2 , B.3| ) according to (|B.1| ) we get 



mHH = {{d-^)a-b)=m (B.4) 



which agrees with ( 3.28 ), and we omitted the superfluous H subscripts ( ^.27|) . 



C. The action in the harmonic ansatz 

Here we derive the form of the action integral (|4.l| ) the action in terms of the harmonic 
fields ( |3.16| ) . Actually the asymptotic region will suffice and there the fields depend only on 
rn- In order to evaluate Ricci's scalar it is useful to have the formula for it in the presence 
of a general fibration 

i 

R = Rx + Yl [e"'^" - 2 di A{Fi) - d,{dFif] -Y,dt dj (dFi ■ dFj) (C.l) 

* id 

where the fibration fields depend only on the x coordinates, Rx, Rvi are the Ricci scalars of 
the spaces X, Yi, Fi = Fi{x), and the Laplacian (A) and grad-squared (d-d) are evaluated 
in the X space. 

Reducing over t, z one has 

= Rd,2 - 2 {dAuf - 2 {dBuf + 2 {BAn ■ OBh) + 2 AAh - 2 ABh (C.2) 

Note that this equation should be (and is) symmetric under the exchange Bh <-> —Ah- 

Since the d — 2 metric is conformal to flat space we may use the expression for the 
Ricci scalar of a conformally transformed metric (see for example p8|) 



■ 2 



ds' = e^'"ds^ 

R = e-^"" [R-2{d-l)Aw-{d-l){d-2) (dw)^] (C.3) 

where d is the dimension of the space and the Laplacian and grad-squared are evaluated 
in the non-tilded metric. In our case w = Ch, d = d — 2 and i? = 0, hence 

Rd-2 = -e-2^« [-2 (d - 3) ACh -{d-3){d- 4) {dCnf] (C.4) 

where the flat space Laplacian can be written as Ad_2 = i9^|^ + {d — 4)/r drH- 
Putting ( |C.2|JC.4[) together one has 



d'^x^dRd = -2dr^^4-3/3oLol7d_3e-^«+^«+('^-^)^« x 

X [{din + '^drH) [-Ah + BH + {d- 3)Ch] + 
+ {d-4)C'H{-A'H + B'H) + 
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where /3o Lq are the periods of the t and z coordinates, the prime denotes a derivative with 
respect to vh and we used A = dj.'jj + {d — 4)/rH dru + (d — 4) C'^ drjj. 
Performing integration by parts we end with 

/ = f dm r'jf^ e-^«+^«+('^-4) ^« K,, XI X' (C.6) 

fed (^N J 

where we use a short-hand notation Xi = [Ah, Bh, Ch] and the matrix Kij is 



K 



-1 -(d-3) 

-1 id-3) 

-(d-3) (d-3) {d-3){d-4) 



(C.7) 



Moreover, the original boundary term in (4.1) is exactly such that it cancels against the 



boundary term generated by the integration by parts so that (p.6| ) is the full action. 
D. Another derivation of Smarr's formula 

Here we give a more direct derivation of the integrated form of the first law (Smarr's 
formula). The free energy is evaluated in terms of the asymptotics, and the derivation 
is completed by combined it with the analogous expression for the mass, and the general 
thermodynamic relation F = E — T S. The latter identity relies on integration by parts to 
relate the horizon with infinity, and thus saves us from performing this step explicitly in 
our derivation. 

The free energy is defined in ( [4.1| ). Since we discuss vacuum black hole solutions we 
have R(i = 0, and so only the second term contributes. It can be rewritten as 

I = -l3F=—^dn[i:^''-^^-'^T.^''~^^] (D.l) 
ottGtv 

where S^"*"^) is the {d - 1) area of the boundary ('^S^'^ ^-^ is the same for the reference 
geometry), and dn is a derivative with respect to orthogonal motion of this boundary. 

Let us compute the free energy in the harmonic ansatz ( p.l6| ) for a boundary at rn = R 
and take the i? — > oo limit. The first term is 

= e-'^^dr,[if3e-^n{Le^n^d-3{rHe''n''-%H=R = 

~ /3 L Orf_3 [(1 - -^) {d - 3) R''-^ + {-aH + bH + {d-3) ch)] = 

= PL Ud-s [{d - 3) i?'^-^ + {bH-aH)]. (D.2) 

The reference geometry is fiat: ds^ = —dt^ + dz^ + dr'^ + dQd-3 with periods Pr = 
P ^-Ah(R) ^ Lfl = Le^^^^^ and r = Re^^'^^^ to match with the boundary. The second 
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term is 



iP e-^«(^)) (L e^«(^)) f7,_3 {d - 3) i?^"^ e^^-^^ 



/3 L J]rf_3 {d - 3) [i?'^-^ + {-an + bH + {d-A) ch)] 



/3LOrf_3 (d-3) R 



d-4: 



(D.3) 



where in the last equality we used the constraint ( |3.18| ). 
Combining ( D.2 , p.3 ) according to ( p.l| ) we get 



PL 



N 



- an] 



F 



(a-b) 



namely 



where we omitted the superfluous H subscripts ( p^ 
Now we may use the relation 

F=m-TS 



(D.4) 
(D.5) 

(D.6) 



where T = k/{2tt) is the temperature, and S = A/{4Gn) is the entropy. This relation 
holds for black hole thermodynamics and is proven using integration by parts. Substituting 



m from ( 3.28 ) we get 



kdGN 



(D.7) 



which coincides with the "numerically adapted form" of Smarr's formula (4.13). 
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